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We study the superconducting correlations induced in graphene when it is placed between two 
superconductors, focusing in particular on the supercurrents supported by the 2D system. For 
this purpose we make use of a formalism placing the emphasis on the many-body aspects of the 
problem, with the aim of investigating the dependence of the critical currents on relevant variables 
like the distance L between the superconducting contacts, the temperature, and the doping level. 
Thus we show that, despite the vanishing density of states at the Fermi level in undoped graphene, 
supercurrents may exist at zero temperature with a natural 1/L 3 dependence at large L. When 
temperature effects are taken into account, the supercurrents are further suppressed beyond the 
thermal length Lt (~ vf/IcsT, in terms of the Fermi velocity vf of graphene), entering a regime 
where the decay is given by a 1/L° dependence. On the other hand, the supercurrents can be 
enhanced upon doping, as the Fermi level is shifted by a chemical potential [i from the charge 
neutrality point. This introduces a new crossover length L* ~ Vf/[j,, at which the effects of the 
finite charge density start being felt, marking the transition from the short-distance 1/L 3 behavior 
to a softer 1/L 2 decay of the supercurrents at large L. It turns out that the decay of the critical 
currents is given in general by a power-law behavior, which can be seen as a consequence of the 
perfect scaling of the Dirac theory applied to the low-energy description of graphene. 

I. INTRODUCTION 

Since the discovery of single atomic layers of carbon in 2004a, this new two-dimensional (2D) material (so-called 
graphene) has attracted a lot of attention^. From the experimental point of view, the 2D carbon sheets have shown 
a number of remarkable electronic properties. Thus, there has been evidence that graphene may have a finite lower 
bound (4e 2 //i) in the conductivity at the charge neutrality points^. Furthermore, an anomalous integer quantum 
Hall effect has been measured in the 2D system with plateaus at odd-integer values of the quantum of conductance^. 
The absence of weak localization effects^ has also pointed at the unconventional effects that impurities and in general 
disorder may produce in the graphene sheet. 

Most of the remarkable transport properties of graphene have to do with its particular band structure at low energies. 
The undoped system has a finite number of Fermi points, placed at the corners of the hexagonal Brillouin zone. Only 
two of such points can be taken as independent, with quasiparticle excitations which have conical dispersion above and 
below the Fermi level^. This explains that the low-energy electronic states of graphene may be accommodated into 
two two-component spinor fields, governed by a Dirac hamiltonian which leads to a dispersion relation e(k) = ±Wj?|k|. 
The electronic system displays hence a relativistic-like invariance at low energies, which is at the origin of the finite 
lower bound in the conductivityi^'iS, the anomalous integer Hall effect^ ' 11 ? 12 , and the absence of backscattering in 
the presence of long-range scatterersii Other exotic effects relying on the Dirac theory have been proposed, like 
the selective transmission of electrons through a n-p junction^ or the specular Andreev reflection at a graphene- 
superconductor interface^. 

Recently, the properties of graphene have been also investigated when the material is placed between superconduct- 
ing contacts. Thus, in the experiment reported in Ref. [T5 , it has been possible to measure supercurrents in graphene 
by attaching wide superconducting electrodes with a spatial separation of w 0.5 /im. In another experiment, reported 
in Ref. [l?], a quite different geometry has been investigated by placing thin electrodes across a large 2D sample, with 
a minimum separation between the tips of « 2.5 /im. In this case, the evidence of the superconducting correlations in 
graphene has been obtained in the form of Andreev reflection peaks in the I-V curves, as well as in the abrupt drop of 
the resistance at a temperature of rj 1 K, below the critical temperature (ss 4 K) of the superconducting electrodes. 
Moreover, supercurrents have been also measured in the experiment reported in Ref. [TR where their development 
may have been favored by the large aspect ratio (~ 10) between the width of the junction and the lead separation (of 
the order of a few hundreds of nanometers). 

It is therefore pertinent to study the way in which the superconducting correlations are induced in graphene when 
it is placed between two superconductors, and how such correlations may depend on the geometry of the experimental 
setup. In this paper we are going to address this issue, focusing in particular on the supercurrents supported by the 
graphene sheet. We will be using a formalism placing the emphasis on the many-body aspects of the problem. This 
will allow us to clarify a number of questions, regarding the dependence of the critical currents on relevant variables 
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like the distance between the superconducting contacts, the temperature, and the doping level of the graphene sample. 
In this respect, our approach can be seen as complementary to that of Ref. 19, where the Josephson effect has been 
studied in terms of Andreev reflection at superconducting contacts, concentrating on junctions with relatively short 
distance between the electrodes. We will be dealing with a framework where the tunneling and propagation of the 
Cooper pairs in graphene play the central role, placing in principle no restriction on the separation that may exist 
between superconducting contacts. 

The content of this paper is distributed as follows. We will set up in section II the formalism needed to describe the 
tunneling and propagation of Cooper pairs in graphene. This will be applied to the computation of the critical currents 
in section III, where we will also discuss the different regimes depending on the interplay between the temperature 
and the distance between superconducting contacts. Section IV will be devoted to extend our analysis to the case of 
finite doping, showing the enhancement experienced then by the supercurrents. Finally, we will summarize our results 
and draw our conclusions in section V. 



II. MODEL OF GRAPHENE JOSEPHSON JUNCTION 

Our purpose is to build a model that incorporates the low-energy properties of electron quasiparticles in graphene 
as well as the tunneling of electrons from graphene to the superconducting electrodes and vice versa. We take into 
account in particular that, below an energy scale of ~ 1 eV, the electron dispersion relation has a conical shape, with 
a dependence of the energy e on momentum k given by e(k) « ifp-lk^. We have to bear in mind that the 2D system 
has actually two independent Fermi points supporting such a conical dispersion, at opposite corners K, —K of the 
hexagonal Brillouin zone. The dynamics of the quasiparticles in graphene can be therefore described in terms of a 
couple of two-component Dirac spinors v]/( ) ; a — 1,2, with a hamiltonia n 20 ! 21 

H = v F J d 2 r *^ a)t (r) cr {a) • 0*W(r) (1) 

where {cr^} are two different suitable sets of Pauli matrices^ (we use units such that K = 1). In the above expression, 
the label of the spinor components is omitted for simplicity, and a sum is taken implicitly over the spin index a as 
well as over the index a running over the two different low-energy valleys of the dispersion. 

The above hamiltonian has to be then complemented with a term accounting for the tunneling of electrons from 
the graphene side to the superconducting electrodes and vice versa. In this respect, we are going to assume that the 
tunneling takes place with equal amplitude for the two sublattices of the graphene honeycomb lattice. This kind of 
junction may be realized in cases where the contacts between graphene and the superconductors preserve the structure 
of the graphene lattice. From a technical point of view, such condition implies that the different spinor components 
and the different low-energy valleys couple with equal amplitude to the superconductors. By denoting the electron 
fields in the respective superconducting electrodes by 'J'si and ^S2j we may write the tunneling hamiltonian for 
contacts along the coordinates x\ = and xi = L as 

H t =Y, t d V *S B)t (^.y)*sj,<r(xi,tf) + h.c. (2) 
i=i,a Jo 

where the parameter t represents the tunneling amplitude. We stress at this point that, while the contacts have a 
width given by W in Eq. @ , the extension of the graphene layer along the transverse y direction is not constrained by 
this parameter in our model. Thus, our description will apply in general to 2D graphene samples, with dimensions in 
both the transverse and the longitudinal direction much larger than the contacts introduced by the superconducting 
electrodes. 

The properties of the superconducting electrodes have to be also incorporated in the model of the Josephson 
junction. For the description of the supercurrents, it will be enough to specify the normal density of states p and the 
order parameter A in the superconducting state. We recall that a supercurrent arises in general from a gradient in the 
phase of the order parameter in a superconductor. In the case of a Josephson junction, the supercurrent is produced 
by a mismatch in the phases Xi an d Xi °f t ne respective order parameters in the superconducting electrodes. The 
Josephson current I s is actually given by the derivative of the free energy with respect to the variable \ — \i ~ Xii 
and it can be therefore expressed as 

7, = 2 e |-A ;B Tlog(Tre-^ T ) (3) 
where T is the temperature and H stands for the full hamiltonian of the model. 
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In order to compute the Josephson current from Eq. ([3]), we will resort to a perturbative expansion in the tunneling 
amplitude t. The structure of the dominant contributions may be however very different depending on the actual 
geometry of the Josephson junction^. In cases where the distance L between the contacts is much smaller than 
the superconducting coherence length £, the supercurrents are built from processes with independent tunneling and 
uncorrelated propagation in graphene of the electrons of a Cooper pair. On the other hand, when L is much larger 
than £, the behavior is governed by the fast tunneling and subsequent propagation of the Cooper pair in graphene, 
as shown schematically in Fig. [TJ This situation corresponds to the case where the time of propagation between the 
contacts is much larger than 1/|A|. Under the assumption of a large |A|, the relevant properties of the superconductors 
may be encoded in the statistical average 

{^>Sj,<r(xj,y,— iTi)^sj,-<r(xj,y,-iT2)) « e lXj p 6(n - r 2 ) (4) 

where the operators are ordered with respect to imaginary time r. 




FIG. 1: Schematic representation of the propagation of Cooper pairs in graphene between two superconductors (SC). 

From inspection of the expansion of the r.h.s. in Eq. ([3]) in powers of the tunneling amplitude, we observe that the 
first nonvanishing contribution to I s appears to fourth order in t, from a statistical average of operators participating 
of the condensates of the two superconductors. The expression of the maximum supercurrent I c (critical current) is 
worked out at that perturbative level in the Appendix, focusing on the regime corresponding to L 3> £. After factoring 
out the relative tunnel conductances at the contacts (given in each case by the dimensionless quantity pt 2 W/vp ), we 
end up with an expression for the behavior of the critical current intrinsic to the 2D graphene layer: 

I? D \T) a 2ev 2 F / d Vl / dy 2 / dr (^(0, y i; 0)*;" a)t (0, yi ; 0)¥< 6) (£, y 2 ; -tr)^- b) (L, y 2 - -«-)) (5) 

Jo Jo Jo 

We observe from |(SJ) that the propagator of the Cooper pairs evaluated over a distance L plays the central role in the 
determination of the supercurrents. We will study in what follows the behavior of this propagator depending on the 
distance L, the temperature, and the doping level. 

III. SUPERCURRENTS AT FINITE TEMPERATURE 

We analyze first the behavior of the supercurrents in graphene when the system is undoped, but is placed at a 
nonvanishing temperature T. The expectation values in the above formulas have to be understood then as statistical 
averages at that finite temperature. The building block for all the calculations is the electron propagator 

G^(r,t) = -i(T*W(0,0)*W+(r,t)> (6) 

This is given in graphene by the propagator for Dirac fermions, in correspondence with the hamiltonian ([1]). In the 
many-body theory at temperature T ^ 0, the imaginary part of that object gets a specific term to account for the 
thermal effects. The full expression of the Dirac propagator becomes in momentum spaced 

G (a) (P,- P ) = yff+fe + «*H + • p)5(-ul + P 2 ) 1 + el l pUkBT (7) 

The Cooper-pair propagator in ([5]) can be computed from the convolution of two Dirac propagators, bearing in mind 
that they correspond to fields at opposite valleys of the graphene dispersion. In doing this operation, we will have to 
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be also consistent with our assumption that the tunneling at the superconducting contacts is the same for the two 
sublattices of the graphene lattice. This means that, when taking the average for the Cooper-pair propagator, we will 

also take a trace in spinor space over the states of the Cooper pairs in sublattice A, given by ^^(k + q) 1 ^ " (— q), 

and in sublattice B, given by \&^|(k + q) 1 ?^ ?( — l)- The Cooper-pair propagator thus defined in momentum space, 
.D(k, lo), can be expressed as 



D(k,u; k ) = iTrJ ^ J A G W (q + k , W9+ ^ )G H U . 



-u q ) (8) 



A nice feature of the diagrammatics of the many-body theory at T ^ is that the terms carrying the dependence 
on temperature do not need to be regularized by means of a high energy cutoff. The contributions at T = 0, however, 
remain finite only when the integrals over the momenta are suitably cut off. In the present model, it is convenient 
to choose a method of regularization of the integrals preserving the relativistic-like invariance of the theory. For this 
purpose, we will adopt an analytic continuation in the number of space-time dimensions 2 ^, that is, carrying out first 
the integrals at general dimension D, and then taking the limit D — > 3. To implement this procedure, we first collect 
the components of the momentum and the frequency to form 3D vectors, q = (vpq, u> q ), k = (vplc, u>k)- Next, we 
may rotate all the 3D vectors to Euclidean space by introducing imaginary frequencies, uJ q = —iuj q . One can easily 
see that the expression of the propagator ([8]) at general dimension D becomes 



L>(k, iiOk)\ T=0 = I dx 
Jo 



1 r AD„ r. 2 o/„2 



d u q 2q 2 ~2k 2 x(l 
(2tt) d ( q 2 + k 2 x(l-x)) 
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1 „ ( D\ 1 „ /_ D 



r 1 - -W )- 7TT7? r 2 - ^ / dxy/k*x(l-x) (9) 



4tt 3 / 2 V 2 / 2 ^ 3/2 V 2 

In the last passage we have made use of standard formulas in dimensional regularization. Quite remarkably, the result 
turns out to be finite in the limit D — > 3. After reverting the rotation back to real frequency, we finally get 

D{^)\ T =o = -iy^"^ (10) 

The part of the Cooper-pair propagator depending on temperature can be computed by using the second term in 
to make the convolution ||5J). For our purposes, we can concentrate on the calculation of the Cooper-pair propagator 
at zero frequency. By adding the result (|10p to the temperature-dependent contribution, we get 

From the results (|10p and (|lip. we can already extract a number of conclusions regarding the behavior of the 
supercurrents in long graphene Josephson junctions. From Eq. ([5]). we can express the critical current for L 3> W as 



4 2D) ( T ) « 2ev 2 F W 2 / — |k| J o (|k|L)^(k,0)e- |k|/fcc (12) 
Jo 2lT 



A short distance cutoff k c has been introduced to regularize the integral over the momentum. This is actually justified 
on physical grounds, since the description of graphene as a continuum in terms of the Dirac theory makes sense at 
distances above the nanometer scale. A sensible choice corresponds to VFk c ~ 1 eV. We will see that, at distances 
such that L > fc^" 1 , the behavior of the critical current is in general not sensitive to the actual value of the cutoff. 

At T = 0, the dependence of the critical current on L can be obtained from the Cooper-pair propagator (| 10[) . 
Actually, we can derive an analytical expression for li (0) by computing the integral in l|12p: 

/>oo 

lf D \<S) ev F W 2 dk |k| 2 J (|k|L)e- |k|/fcc (13) 

= e^4ffi£=a (14) 

V(fc c 2 £ 2 + l) 5 

(15) 

From this result we check that, as expected, the behavior of the critical current is not affected by the cutoff k c in the 
limit of large L. In this regime we find 

l( 2 °\0)~ev F W 2 ± (16) 
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FIG. 2: Logarithmic plot of the critical current ll? D ^ (in units of 10 2 evpk c ~ 1.2 fiA) as a function of the distance L, taking 
W = 10 2 /k c (= 50 nm). The three curves correpond, from top to bottom, to different values of the temperature T = 2 K, 4 K, 
and 8 K. 



The strong power-law decay shown by (|16p can be understood actually as a reflection of the linear dependence on 
momentum of the quasiparticle energy, which dictates in turn the behavior of the Cooper-pair propagator (|10[1 25 . We 
reach anyhow the interesting conclusion that, while graphene has a vanishing density of states at the Dirac point, it 
may still support a nonvanishing supercurrent when the Fermi level is at that charge neutrality point. 

The inspection of the full propagator (TlTj) also reveals that the scaling is drastically modified when ksT 3> vj?|k|. 
Actually, we can distinguish between a high-temperature and a low-temperature regime of the Cooper-pair propagator, 
with quite different behaviors: 

£>(k,0) a -^-|k| if k B T<^v F \k\ (17) 

OVp 

log(2) 1 |k| 2 

~ ~ k B T ~ TTTT^? lf k B T^>v F \\i\ (18) 

TTV F l07T Kb-L 

The existence of this crossover in the momentum gives rise to an abrupt decay of the supercurrent beyond the thermal 
length Lt — vp/ksT. This is illustrated in Fig. [21 where the critical current I^ D \t) is represented as a function of 
the distance L at different temperatures. We observe for instance that, for a temperature of the order of T ~ 1 K, the 
scale of the crossover in L is of the order of a few microns, in agreement with the expression of the thermal length. 

From a physical point of view, it becomes clear that the Cooper pairs do not feel the thermal effects during their 
propagation when L is shorter than the scale given by Lt, while they are increasingly disrupted at distances larger 
than the thermal length. At short distances such that L <C vp/ksT, the decay of the critical current represented in 
Fig. [2]follows a 1/L 3 power-law, in agreement with the above analysis at T = 0. However, beyond the crossover clearly 
identified in the three curves, we see that a different power-law behavior opens up at long distance L 3> vp/ksT. 
This regime can be analyzed by considering that, when T is very large, the second term in the approximation (|18p 
dictates the long-distance decay of the critical current. In this case we can compute again analytically the integral in 

I? D \T) ~ -ev 2 F W 2 — J dk |k| 3 J (|k|L)e-l k l/^ (19) 
= ev F W^f k ° L2 -V (20) 

k B T ^UTTy 

The leading contribution to the critical current becomes then for L VF/ksT 

* a ><n~~»"'i%srh < 21 ' 

The existence of this stronger power-law decay is manifest in the results of the numerical computation of the critical 
current represented in Fig. [21 as it can be checked that the rightmost part of the lower curves in the plot corresponds 
with great accuracy to a power-law behavior with the exponent given by Eq. (|2"Tj) . 
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FIG. 3: Plot of the critical current li 2D) (T) (in units of 10 -2 evFk c ~ 1.2 fiA) as a function of the temperature, for W = W 2 /k c 
(= 50 nm) and a spatial separation between superconducting contacts L = 0.5 fim (a), 1.5 /im (b), and 2.5 /im (c). 




FIG. 4: Combined plot of 7 C (T)// c (0) represented as a function of the scaled variable T/T* (with T* = vf/2L), where 
it is seen the collapse of the three curves corresponding to values of the distance L — 0.5 (im (full line), 1.5 /jm (dotted line), 
and 2.5 /J,m (dashed line). 



In order to establish a comparison with experimental results, the relevant behavior is given by the critical current 
represented as as function of the temperature at fixed length L. The existence of a thermal length has a reflection 
here in the form of a crossover temperature T* , which marks the strong decay of the critical current for T > T* . We 
have plotted in Fig. [3] the critical current Iq 2D \t), computed from Eq. (fT2")l , at different values of L between 0.5 /im 
and 2.5 fim. The shapes of the curves in the figure are quite similar, and it can be checked that they can be collapsed 
into a single universal curve after rescaling the temperature by T* oc vp/ksL, as shown in Fig. [4] This is consistent 
with the expression of the critical current in Eq. (fT2|) , where it it seen that the effect of a variation of the length L on 

Ic 2D \t) / Ic 2D \o) can be compensated by a suitable change in the scale of T, in the regime where the critical current 
is not sensitive to the precise value of k c . 

We observe that the behavior of the critical current is in all cases quite stable for T <C T* and that there is even an 
upturn before the abrupt drop at the crossover temperature. These features have been also found in the theoretical 
investigation of the supercurrents in one-dimensional (ID) electron systems^ and in carbon nanotubes^.. The shape 
of the critical currents obtained there is qualitatively similar to that of the curves in Fig. [3J A major difference is 
however that the decay of the supercurrents in the carbon nanotubes is given by a 1/L dependence in the ballistic 
regime, instead of the much stronger power-law decay (| 16|) in graphene. 

It is worth mentioning at this point the experiment reported in Ref. UtI in which the properties of a graphene 
Josephson junction have been measured in the regime of large distance between superconducting electrodes. In the 
experimental setup described there, the minimum distance between superconducting contacts can be estimated as 
« 2.5 ^m. While no supercurrent was observed below the critical temperature T c of the electrodes (« 4 K), a signature 
of the proximity effect was obtained in the measurements of the resistance as a function of temperature, in the form 
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of a sharp decrease at T ~ 1 K. Quite remarkably, this value of T is in good correspondence with the crossover 
temperature that we find in our model for a distance L — 2.5 /im, as can be seen from Fig. Efc). It is therefore 
likely that the sharp decrease measured in the resistance has its origin in the same suppression of the thermal effects 
that enhances the supercurrents at T < T*. We also notice that the prediction from our model is that the critical 
currents for such a large value of L should be well below the scale of 1 nA. This may explain the failure to establish a 
supercurrent in the experiment of Ref. ITtI and it may also anticipate better perspectives in experiments with suitably 
short graphene junctions. 

IV. SUPERCURRENTS AT FINITE DOPING 

We have seen that the origin of the relative smallness of the critical currents in undoped graphene lies in the 
vanishing density of states at the Dirac point. Therefore, a straightforward way to enhance the supercurrents may 
simply consist in shifting the Fermi level away from the charge neutrality point, as shown in Fig. [5l In practice, this 
can be achieved by doping the graphene sheet. In our theoretical framework, we will assume that this effect can be 
accounted for by means of a finite chemical potential fi . Thus, the hamiltonian for the graphene part of the junction 
will now read: 

H = [ d 2 r *i a)t (r) (v F (T^ ■ d - y) 1^ Q) (r) (22) 





FIG. 5: Schematic representation of the two independent Dirac valleys at the corners of the hexagonal Brillouin zone, showing 
the regions of occupied (dark) and unoccupied (white) energy levels in doped graphene. 

Working at fi ^ leads to significant modifications in the propagator of the Dirac fermions and in the Cooper-pair 
propagator. The Dirac propagator corresponding to the hamiltonian (|22[) turns out to be (for [i > 0)21 

G( Q )(k,c) = (u + v F crW-k){- 2 Lf— + m 5{u ~ "/^W - v F \k\)] (23) 

w — v F K^ + ie vf\k\ 

As shown in the Appendix, the representation (|23p is nothing but a compact form of expressing the propagation of 
quasiparticles with Wp|k| > fi and quasiholes with ±u^|k| < /i, in the particular case of conical dispersion. 

The propagator (|23|) is very convenient to carry out calculations in the many-body theory and, in particular, it 
allows us to compute the dependence on /i of the Cooper-pair propagator as a correction to the expression (I10p at 
/j = 0. In this procedure, we observe that the second term in the r.h.s. of (f2"3"]l does not introduce any integrals 
requiring regularization in the diagrammatics of the Dirac theory. By computing then the Cooper-pair propagator 
according to Eq. flSJ, we obtain^ 



D(k,0) = -tt^tM if v F \k\ < 2y 

iTTVp 

- -i' k ' + 4^ |k|arcsin (^ki) - ^ if ^ |k| > ^ (24) 

At large values of Vp\k.\ ^> n, we recover from (pM)) the linear dependence on the momentum that is characteristic of 
the Cooper-pair propagator in the undoped system. However, the chemical potential introduces a clear deviation from 
that behavior at small |k|, which has significant consequences in the decay of the supercurrent at long distances. This 
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FIG. 6: Plot of the zero-temperature critical current 7^ (0) as a function of the distance L, for W = 10 2 /fc c (= 50 nm) and 
three different values of the chemical potential [i = 1 meV (a), 5 meV (b), and 10 meV (c). The dashed straight lines in Fig. 
3(a) are drawn as a reference to the power-law dependences 1/L 3 and 1/L 2 . 



is illustrated in Fig. [5J where it can be appreciated the existence in general of a crossover length scale L* mediating 
the transition towards a softer power-law decay. 

According to (pM]) , we can express the critical current 1^?°^ (0) at finite chemical potential in the form 



j(^) ( o)=ir ) (o)+4 2D) (o) 



with 



4l D \0) = -ieW% r^|k|J (|k|L )e -l k ^ 

7T Jq ZTT 



4 2D) (o) 



2tt 



ev F W 2 



dk 

2^ 



Ikl 2 arccos 



2fi 



v F \k\ 



J (|k|L) e -l k l/^ 



The first contribution to (|25|) is not relevant, since we have 



4 D) (o) 



-eW 2 (i 



(25) 

(26) 
(27) 

(28) 



which is smaller than the estimate (fTB| at fi = by a factor n/v F k c . The second contribution may change however 
the behavior of the critical current at large L, as the integrand is not analytic at |k| = 2/i/vp. The expression for 
I& D \o) is actually finite in the limit k c — > oo, and we obtain 

,3 



(2D) 



(0) 



-ev F W 2 ^- 



dx 



arccos | — ) Jq{{2^lL/vf)x) 



for [ihjvF ^S> 1 



(29) 
(30) 



We have to stress anyhow that I^ D> {Q) has oscillations has a function of L, arising from the own behavior of the 
Bessel function J . The power-law decay ([50)) applies then to the envelope of the maxima of the critical current, as it 
is illustrated in Fig. [6^a). There it can be appreciated the crossover from the 1/L 3 behavior to the oscillatory regime 
with softer power-law decay. From the numerical results represented in the figure, it can be checked that the 1/L 2 
behavior is followed with great accuracy at larges values of L (compared to vp/fi). 

From a practical point of view, the most important result that we obtain is the significant enhancement of the 
critical currents at moderate values of the chemical potential. This is clearly observed in the plots of Fig. [6l where 
the crossover to the 1 /L 2 decay is always found at a length scale consistent with the theoretical estimate L* ~ vf /m- 
For a chemical potential /i ss 10 meV, for instance, that scale is » 50 nm. The critical currents can be then enhanced 
to values above the nanoampere scale for spatial separation between superconducting contacts L > 500 nm (assuming 
thin electrodes as in our case with W ~ 50 nm). This should open good perspectives to establish supercurrents in 
graphene Josephson junctions by suitable doping of the samples. 
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V. CONCLUSION 

In this paper we have adopted a framework suited to address the many-body properties of graphene Josephson 
junctions. We have described the development of the supercurrents through the tunneling and propagation of Cooper 
pairs in the graphene part of the junction, with the aim of investigating the dependence of the critical currents on 
relevant variables like the distance between the superconducting contacts, the temperature, and the doping level. We 
have been able then to characterize different regimes in the behavior of the supercurrents, depending on the relation 
between those variables. 

The supercurrents have a natural tendency to decay in the graphene part of the Josephson junction, following in 
general a power-law behavior with respect to the distance L between the superconducting contacts. Such a power- 
law decay is particularly strong in undoped graphene, given the vanishing density of states at the charge neutrality 
point. We have shown that the critical currents display then at zero temperature a 1/L 3 dependence on the distance 
L. When temperature effects are taken into account, there is always a finite thermal length Lt (of the order of 
~ vp/ksT) beyond which the supercurrents are further suppressed, due to the disruption of the Cooper pairs by 
many-body effects. When this takes place, the supercurrents enter a regime where the natural decay is given by a 
1/L 5 dependence. 

On the other hand, many-body effects can be also used in our benefit to enhance the critical currents, in this case 
by shifting the Fermi level away from the charge neutrality point. This can be achieved in our framework by means 
of a chemical potential /.i ^ 0. Inducing in this way a finite density of states at the Fermi level, we have seen that 
the critical currents are enhanced beyond a new crossover length L* ~ Vf//J- This is actually the scale at which the 
effects of the finite charge density start being felt, marking the transition from the previously discussed 1/L 3 behavior 
to a softer 1/L 2 decay of the supercurrents at long distances. 

At this point, it is interesting to note that the 1/L 3 decay at zero temperature in undoped graphene is similar to the 
behavior found in the investigation of mesoscopic junctions made of a diffusive metal 2 -. In this case, the product of 
the critical current times the normal resistance of the metal is proportional to the Thouless energy, which depends on 
length L as 1/L 2 . This implies consequently a 1/L 3 decay of the critical current, which we have seen is characteristic 
of graphene under conditions of ballistic transport. The reminiscence of some of the properties of clean graphene with 
respect to the behavior of a disordered normal metal has been remarked in several other instance s] 9 ' 19 . We have to 
point out, however, that this resemblance does not go farther in our case, regarding other regimes of the graphene 
Josephson junction. In particular, we have seen that the critical current does not follow an exponential decay at 
distances larger than the thermal length. The decay of the critical current is always given in graphene by a power-law, 
which can be seen as a consequence of the perfect scaling of the low-energy Dirac theory. 

We have also to stress that our results refer to Josephson junctions with graphene layers which have large dimensions 
in both the longitudinal direction along the junction and the transverse direction. This condition comes from our 
consideration of a system which is truly 2D, where in particular the size in the direction transverse to the junction 
is not constrained by the width of the superconducting contacts. In this regard, the situation is quite different to 
the case of long but narrow junctions, where the small transverse dimension may lead to the quantization of the 
transverse component of the momentum. In such circumstances, the behavior of the system may be rather dictated 
by a ID propagation of the Cooper pairs, which is known to lead to a 1/L decay of the supercurrents in the ballistic 
regim o 22 ! 26 . 

Anyhow, the great advantage of the graphene Josephson junctions is that the interaction effects have little signifi- 
cance at the temperatures required to measure the supercurrents. In the long ID junctions made of carbon nanotubes, 
for instance, it is known that the Coulomb interaction may induce a strong power-law suppression of the density of 
states, with the consequent reflection in the decay of the supercurrent^. In 2D graphene, however, the electron system 
has the tendency to become less correlated at low energies, with a strong renormalization of the Coulomb interaction 
that makes it practically irrelevant at the temperature scale of 1 K 24 i 29 . 

In conclusion, our results highlight the role of the different parameters conforming the geometry of graphene 
Josephson junctions in the determination of the critical currents. We have seen that the interplay with variables like 
the temperature and the doping level is what establishes the different regimes of a junction. This information may be 
useful in the design of experiments, for the purpose of enhancing the magnitude of the critical currents in real devices. 

Acknowledgments 

We thank H. Bouchiat and F. Guinea for useful comments. The financial support of the Ministerio de Education 
y Ciencia (Spain) through grant FIS2005-05478-C02-02 is gratefully acknowledged. E.P. is also financially supported 
by CNISM (Italy). 



10 



Appendix 

A. Lowest-order contribution to the critical current 

The Josephson current I s can be computed in a perturbative framework by expanding the free energy in Eq. ([3} in 
powers of the tunneling amplitude t. The first nonvanishing contribution is found to fourth order in this expansion, 
as the statistical average of operators leads then to the appearance of the condensates of the two superconductors in 
the junction. We have actually 

q 4 .\/k B T rW 

I, ~ 2e— fc B r« 4 [] | dnj dy l (^ slA (0,y 1 ;-iT 1 )^ sld {0,y 2 ;-iT 2 )) 

x (*<»»((), yi; -in)*i~ o)t (0, y 2 ; -ir 2 )^ (L, y 3 ; -in)9^ b) (L, y 4 ; -ir 4 )) 

x (^kt^' y3 ' - ir 3)*L,i( L ' - iT ^>) ( 31 ) 

We can apply to Eq. (|31[) the approximations pertinent to the regime we want to study in the paper. Focusing on 
the case of a large junction where the distance L is much larger than the superconducting coherence length £, the use 
of Eq. {1} and translational invariance in the variable r leads to a maximum supercurrent I c (critical current) 

2 pW fl/k B T 

I C (T) « 2ep 2 t i W 2 H / % / ^(^(O.ifcjO)^"""^^^^ (32) 

i=1 Jo Jo 

At this point, it becomes convenient to factor out the relative tunnel conductances at the contacts, which are each 
given by the dimensionlcss quantity pt 2 W/vF- We concentrate then on the behavior of the critical current intrinsic 
to the 2D graphene system, represented by the expression 

lf D \T) » 2evl / dyj dy 2 dr <^ o)t (0, Vl ; 0)*;" a)t (0, m 0)* f\ L , y 2 ; -ir)¥<- 6) (L, y 2 ; -ir))(33) 

Jo Jo Jo 

This last equation highlights the connection between the critical current and the propagator of the Cooper pairs, 
which plays a central role in the discussion of Sections III and IV in the paper. 

B. Dirac propagator at fj, ^ 

In the many-body theory of Dirac fermions, it is usual to write the Dirac propagator at finite charge density in the 
form (assuming fi > 0) 

GW(k, W ) = (c + W a) -k)[- 1-—- + wr — ^— ~ gff U fa - v F \k\)} (34) 

lo z — v F \s. z + ie VF\k\ 

If we specialize the expression (jMJ) to modes such that the eigenvalue e(k) of the matrix vpcr^ ■ k is positive, we 
get 

WaVl \ 1 . U! + V F \k\ 

G ( ) (k,w) = prr -«r ' ' 5 (u - v F \k\) 

' (35) 



for uj?|k| > /i, and 



w — vp\k\ + ie 



G«(M) = _J_+^i^±^lff(a,-»,,|k|) 

e(k)=« F |k| W-U_F|k| 2v F |k| 



= TUT" T ( 36 ) 

for Vp\k\ < (j,. We observe that, in either case, the expression of G^°'(k, oj) agrees with the conventional propagation 
of a fermion, with the correct ±ie prescription depending on whether it corresponds to a quasiparticle or a quasihole 
excitation. 
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On the other hand, we always have in the case of negative eigenvalue e(k) 

1 lu — vf\\c\ 



G (a) (k,c^) 



e(k)=-v F \k\ W + W F |k| 
1 



2v F \k\ 



5(u; + v F \k\) 



which is also in agreement with the expected propagation for a quasihole in the valence band of graphene. 
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